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Abstract. The presence of convective motions in the atmospheres of metal-poor halo stars leads 
to systematic asymmetries of the emergent spectral line profiles. Since such line asymmetries are 
very small, they can be safely ignored for standard spectroscopic abundance analysis. However, 
when it comes to the determination of the 6 Li/ 7 Li isotopic ratio, q(Li)=n( 6 Li)/n( 7 Li), the 
intrinsic asymmetry of the 7 Li line must be taken into account, because its signature is essentially 
indistinguishable from the presence of a weak 6 Li blend in the red wing of the 7 Li line. In 
this contribution we quantity the error of the inferred 6 Li/ 7 Li isotopic ratio that arises if the 
convective line asymmetry is ignored in the fitting of the A 6707 A lithium blend. Our conclusion 
is that 6 Li/ 7 Li ratios derived by Asplund et al. (2006), using symmetric line profiles, must be 
reduced by typically Ag(Li) ~ 0.015. This diminishes the number of certain 6 Li detections from 
9 to 4 stars or less, casting some doubt on the existence of a 6 Li plateau. 
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1. Introduction 

The spectroscopic signature of the presence of 6 Li in the atmospheres of metal-poor 
halo stars is a subtle extra depression in the red wing of the 7 Li doublet, which can 
only be detected in spectra of the highest quality. Based on high-resolution, high signal- 
to-noise VLT/UVES spectra of 24 bright metal-poor stars, Asplund et al. (2006) report 
the detection of 6 Li in nine of these objects. The average 6 Li/ 7 Li isotopic ratio in the 
nine stars in which 6 Li has been detected is <z(Li) w 0.04 and is very similar in each 
of these stars, defining a 6 Li plateau at approximately logn( 6 Li) = 0.85 (on the scale 
logn(H) = 12). A convincing theoretical explanation of this new 6 Li plateau turned 
out to be problematic. Even when the depletion of the 6 Li isotope during the pre-main- 
sequence phase would be ignored, the high abundances of 6 Li at the lowest metallicities 
cannot be explained by current models of galactic cosmic-ray production (for a concise 
review see e.g. Christlieb 2008, and references therein). 

A possible solution of the so-called 'second Lithium problem' was suggested by Cayrel 
et al. (2007), who point out that the intrinsic line asymmetry caused by convection in the 
photospheres of cool stars is almost indistinguishable from the asymmetry produced by 
a weak 6 Li blend on a presumed symmetric 7 Li profile. As a consequence, the derived 6 Li 
abundance should be significantly reduced when the intrinsic line asymmetry in properly 
taken into account. Using 3D non-LTE line formation calculations based on 3D hydro- 
dynamical model atmospheres computed with the C0 5 BOLD code (Freytag et al. 2002, 
Wedemeyer et al. 2004, see also http://www.astro.uu.se/~bf/co5boldjnain.html), 
we quantify the theoretical effect of the convection-induced line asymmetry on the re- 
sulting 6 Li abundance as a function of effective temperature, gravity, and metallicity, for 
a parameter range that covers the stars of the Asplund et al. (2006) sample. 
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2. 3D hydrodynamical simulations and spectrum synthesis 

The hydrodynamical atmospheres used in the present study are part of the CIFIST 3D 



model atmosphere grid, as described by Ludwig et al. (2009) They have been obtained 
from realistic numerical simulations with the CO b BOLD code which solves the time- 
dependent equations of compressible hydrodynamics in a constant gravity field together 
with the equations of non-local, frequency-dependent radiative transfer in a Cartesian 
box representative of a volume located at the stellar surface. The computational domain 
is periodic in x and y direction, has open top and bottom boundaries, and is resolved 
by typically 140x140x150 grid points. The vertical optical depth of the box varies from 
logTRoss ~ —8 (top) to logTR OSS w +8 (bottom). The selected models cover the stellar 
parameter range 5900 K < T cff < 6500 K, 3.5 < log 5 < 4.5, -3.0 < [Fe/H] < -1.0. 

Each of the selected models is represented by about 20 snapshots chosen from the full 
time sequence of the corresponding simulation. All these representative snapshots are 
processed by the non-LTE code NLTE3D that solves the statistical equilibrium equa- 
tions for a 17 level lithium atom with 34 line transitions, fully taking into account 
the 3D thermal structure of the respective model atmosphere. The photo-ionizing ra- 
diation field is computed at 704 frequency points between A 925 and 32 407 A, using 



the opacity distribution functions of Castelli & Kurucz (2004) to allow for metallicity 



dependent line-blanketing, including the HI-H+ and HI-HI quasi-molecular absorption 
near A 1400 and 1600 A, respectively. Collisional ionization by neutral hydrogen via 
the charge transfer reaction H(ls) + Li(n£) <-> Li + (ls 2 ) + H~ is treated according to 



Barklem, Belyaev & Asplund (2003) More details are given in Sbordone et al. (2009) 

Finally, 3D non-LTE synthetic line profiles of the Li I A 6707 A feature are computed 
with the line formation code Linfor3D (http://www.aip.de/^mst/linfor3D_main.html), 
using the departure coefficients bi = rii(NLTE)/ni(LTE) provided by NLTE3D for each 
level i of the lithium model atom as a function of geometrical position within the 3D 
model atmospheres. As demonstrated in Fig.[TJ 3D non-LTE effects are very important 
for the metal-poor dwarfs considered here: not only is the 3D LTE equivalent width too 
large by more than a factor 2, but also is the half- width of the 3D LTE line profile too 
narrow by about 10%. Moreover, the lithium lines are significantly less asymmetric if the 
non-LTE effects are taken into account. 



3. Method and Results 

As outlined above, the 6 Li abundance is necessarily overestimated if one ignores the 
intrinsic asymmetry of the 7 Li line profile. To quantify this error theoretically, we rely 
only on synthetic spectra. The idea is to represent the observation by the synthetic 
3D non-LTE line profile of the 7 Li line blend. This 3D flux profile is computed with 
zero 6 Li content. Except for an optional rotational broadening, the only source of non- 
thermal line broadening is the 3D hydrodynamical velocity field, which also gives rise to 
a convective blue-shift and an intrinsic line asymmetry. Next we compute a small grid of 
ID LTE synthetic line profiles of the full 6 Li/ 7 Li blend from a so-called ID LHD model, 
a ID mixing-length model atmosphere that has the same stellar parameters and uses the 
same microphysics and radiative transfer scheme as the corresponding 3D model. The 
parameters of the grid are the total 6 Li+ 7 Li abundance, A(Li), and the 6 Li/ 7 Li isotopic 
ratio, <?(Li). Microturbulcnce is fixed at £ m ; c = 1.5 km/s, v sini is identical to the value 
used in the 3D spectrum synthesis (we tried and 2 km/s). Now the ID line profiles from 
the grid are used to fit the 3D profile. Four parameters are varied independently to find 
the best fit (minimum \ 2 ): in addition to A(Li) and g(Li), which control line strength 
and line asymmetry, respectively, we also allow for an extra line broadening characterized 
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Figure 1. Comparison of 3D LTE (dashed) and 3D non-LTE (solid) line profile (left) and line 
bisector (right) of a single 7 Li component, computed for a typical metal-poor turn-off halo star 
(T eff = 6300 K, log g = 4.0, [Fe/H] = -2). Non-LTE effects strongly reduce the equivalent 
width of the line (W: 35.5 — * 15.6 mA) while they increase the half width of the line profile 
[FWHM: 7.7 — > 8.5 km/s); the line asymmetry is diminished in non-LTE (velocity span of 
bisector 5v : 0.31 -> 0.19 km/s). 



Table 1. Fitting three observed Li I A 6707 A spectra with ID LTE and 3D non-LTE synthetic 
line profiles, respectively. Columns (7)-(10) show the results for v sin i — 0.0 / 2.0 km/s. 
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HD 74000 


6203 


4.03 


-2.05 


600 


3D NLTE 


2.25 / 2.25 


-1.1 / -1.1 


3.1 / 2.1 


0.64 / 0.64 












ID LTE 


2.23 / 2.23 


0.6 / 0.6 


5.9 / 5.4 


0.42 / 0.42 


G271-162 


6230 


3.93 


-2.30 


550 


3D NLTE 


2.30 / 2.30 


0.6 / 0.5 


3.6 / 2.9 


0.04 / 0.05 












ID LTE 


2.27 / 2.27 


2.2 / 2.2 


6.2 / 5.7 


-0.17 / -0.17 


HD 84937 


6310 


4.10 


-2.40 


630 


3D NLTE 


2.21 / 2.20 


4.0 / 4.2 


3.7 / 2.7 


0.08 / 0.07 












ID LTE 


2.18 / 2.18 


6.3 / 6.0 


6.1 / 5.6 


-0.17 / -0.14 



Notes: 11 T off /logs/ [Fe/H] = 6280K/4.0/-2; 2) log [n( 6 Li) + n( 7 Li)] - logn(H) + 12; 3> Gaussian kernel 

by FWHM of the Gaussian kernel, and a global line shift, Av. The value Q*(Li) of the 
best fit is then identified with the correction Aq(Li) that has to be subtracted from 
the 6 Li/ 7 Li isotopic ratio determined from the ID LTE analysis in order to correct for 
the bias introduced by the intrinsic line asymmetry: Ag(Li) = g*(Li), and q^ 3D ^(Li) = 
g( 1 - D )(Li) - Aq(Li). The procedure takes saturation effects properly into account. 

We have determined A q(Li) in the relevant range of stellar parameters according 
to the method outlined above. The results are displayed in Fig. [5] for [Fe/H]=— 1 and 
—2. At given metallicity, the corrections are largest for low gravity and high effective 
temperature, increasing towards higher metallicity. We note that A <?(Li) is essentially 
insensitive to the choice of v sini. The downward correction of the 6 Li/ Li isotopic ratio 



is typically in the range 0.01 < Ag(Li) < 0.02 for the stars of the Asplund et al. (2006) 
sample (see Fig. [5]). After subtracting for each of these stars the individual Ag(Li), ac- 
cording to T e g, \ogg, and [Fe/H], the mean 6 Li/ 7 Li isotopic ratio of the sample is re- 
duced from 0.0212 to 0.0059, as illustrated in Fig. 3, If we keep the error bars given by 
Asplund et al. (2006) . the number of stars with a 6 Li detection above the 2tr level de- 



creases from 9 to 4. One of them, HD 106038, survives only because of its particularly 
small error bar of ±0.006, another one, CD-30 18140, just barely fulfills the 2a crite- 
rion. The remaining two stars are G020-024, which shows the clearest evidence for the 
presence of 6 Li (q(Li) = 0.052 ± 0.017), and HD 102200 with a somewhat weaker 6 Li 
signal (q(Li) = 0.033 ± 0.013). The spectra of these stars should be reanalyzed with 3D 
non-LTE line profiles. 

As a consistency check, we have also fitted a few observed Li I A 6707 A spectra with 
ID LTE and 3D non-LTE synthetic line profiles, respectively. The fitting parameters are 
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again A(Li), <?(Li), FWHM, and Aw. As expected, the 3D analysis yields lower q(Li) by 
roughly -0.02. Details are compiled in TableQ] HD 74000 and G271-162 are considered 
non-detections, while HD 84937 remains a clear 6 Li detection with q( 3£l )(Li) w 0.04. 

4. Conclusions 

The present study indicates that only 2 or at most 4 out of the 24 stars of the 



Asplund et al. (2006) sample remain significant 6 Li detections when subjected to a 3D non- 
LTE analysis, suggesting that the presence of 6 Li in the atmospheres of galactic halo stars 
is rather the exception than the rule. This would imply that it is no longer necessary to 
look for a global mechanism accounting for a 6 Li enrichment of the galactic halo, but 
that it is sufficient to explain only a few exceptional cases, which is probably much easier. 
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Figure 2. Contours of Aq(Li) in the T e s - \ogg plane for metallicities 
[Fe/H] = — 1 (right). White symbols mark the positions of the stars from the 
sample with -2.5 < [Fe/H] < -1.5 (left), and -1.5 < [Fe/H] < -0.5 (rig 



Fe/H] = -2 (left) and 
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Figure 3. 6 Li/ 7 Li isotopic ratio, and ±la error bars, as a function of effective temperature as 



derived by Asplund et al. (2006) before (left) and after (right) subtraction of Ag(Li) to correct 
for the bias due to the intrinsic line asymmetry. Filled circles denote 6 Li detections above the 
2cr level, open circles denote non-detections. 



